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This talk summarizes the status of the calculations of Bk, Br, Bg, and B„, done in collaboration with T. 
Bhattacharya, G. Kilcup, and S. Sharpe. Results for staggered, Wilson, and Clover fermions are presented. 



1. INTRODUCTION 

Reliable estimates of the matrix elements of 
4-fermion operators between hadronic states are 
essential in order to quantify strong interaction 
corrections to weak processes. Here we report 
on calculations of mixing elements of (i) As = 2 
operators between K° and K° states (Bk) that 
arise in the calculation of the CP violation pa- 
rameter e. The value of Bk is an essential in- 
put in pinning down the Wolfenstein parameters 
p and T) in the CKM matrix, (ii) The strong and 
electromagnetic penguin operators needed to pre- 
dict e'/e (B 6 and B$). In particular, we consider 
Bg, which is phenomenologically important since 
a smaller value means a larger e'/ e. (iii) The S+P 
operators needed in the study of the lifetime dif- 
ferences of B mesons (B s ). 

2. STAGGERED RESULTS 

Staggered fermions are the method of choice for 
calculating kaon matrix elements as they respect 
the continuum chiral Ward identities. Our results 
for Bk,By, and Bg have recently been given in 
(lj. They are based on the same numerical data 
first presented at LATTICE 93 by Sharpe @. As 
explained in GJ, even though the statistical qual- 
ity of the data is meager by present standards, the 
largest source of error in the quenched staggered 
theory is a systematic one - the dependence of Bk 
on the lattice operator. The second new feature 
of the analysis, compared to Q, is a better under- 
standing of the matching between lattice and con- 
tinuum operators using tadpole improved 1-loop 
perturbation theory, i.e. the horizontal matching 



explained in ||. 

The data for both "smeared" (/3 = 6.0,6.2) 
and "unsmeared" (j3 = 6.0, 6.2, 6.4) operators are 
shown in Fig. [lj (See [|l| for definition of these 
operators.) The extrapolation to a — is done 
assuming that only the leading correction 0(a 2 ) 
contributes (the absence of the linear 0(a) term is 
expected theoretically || , and has been confirmed 
numerically Q). The results are (for q* = 1 /a) 

Bk(NDR,2 GeV) = 0.63(2) unsmeared 
B K (NDR,2GeV) = 0.60(2) smeared (1) 



For the final value we take the mean 
B K (NDR, 2 GeV) = 0.62 ± 0.02 ± 0.02 



(2) 
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where the second error covers the spread due to 
the operator dependence. 

The difference between the central values for 
smeared and unsmeared operators is ~ 0.025. 
While our data is certainly not good enough to 
argue that this difference is significant, a differ- 
ence of similar size has been reported by JLQCD 
Q . This difference is an artifact of keeping only 
an a 2 correction term in the a = extrapolation 
as shown by the following argument. Consider 
two discretizations, 0\ and 2 , of any 4-fermion 
operator. Let the typical lattice momenta q* as- 
sociated with their lattice measurement be K\/a 
and K 2 1 a respectively. Then, using the "horizon- 
tal" matching to the continuum scheme denned in 
|J, the results at /1 = K\/a are related as 

0x(m) = 2 ^)(a(K 2 /a)/a(K 1 /a)r^/ 2 ^ 

+ {1+Xa 2 (q*)+Ya 2 + ...} . (3) 

Since the factor a(K 2 / a) / a(K\ / a) — *• 1 as a — * 1, 
the two operators should give the same result in 
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Figure 1. The data for B K { NDR, 2 GeV) as a 
function of lattice spacing a 2 , along with a lin- 
ear extrapolation to a = 0, for the smeared and 
unsmeared operators. 



mixing contribution due to the P ® P term can 
be even larger than the tree level result. In such 
a case the 1-loop determination of the matching 
Z's is inadequate, and no results for B7 or B$ 
can be extracted. This is true for the unsmeared 
Landau gauge operators that we have used and 
probably also for the gauge invariant operators 
used by JLQCD Q as they have the same Zp. 
On the other hand, the 1-loop perturbative value 
for Zp for the smeared operators is much smaller, 
and consequently results are independent of q* to 
within 10% as shown in Table [I]. (The failure 
of 1-loop Z's for the unsmeared operators shows 
up as a large dependence of B-parameters on q*\ 
for example the results even change sign between 
q* = ir/a and 1/a @.) 

The shaky part of this analysis is that the data 
at two values of a are extrapolated using just the 
lowest order (a 2 ) correction. Since these correc- 
tion are large, further checks of these first results 
are needed. 



the continuum limit provided the extrapolation is 
done including factors of both 0(a) and O(a). 

To convert the quenched result in Eq. || to the 
renormalization group invariant quantity B, one 
can proceed in two ways: use the rif = or the 
nf = 4 values for a s , /?o, 01, 71 • The quenched 
result, using a(2 GeV) = 0.19, is 

B K = 0.86 ±0.03 ±0.03 (4) 



Interestingly, using n/ = 4 values for /3o, /?i , 71 

and o^t(2 GeV) = 0.3 also gives the same value. 
However, there is an uncertainty of ~ 0.05 in 
such a conversion. What we really predict is the 
quenched result Eq. ^|. 

2.1. Electromagnetic penguins 

The quantity that enters in the standard model 
calculation of e'/e is B$ evaluated at say /j, — 2 
GeV H . To get this from lattice simulations one 
needs to calculate the matrix elements of both O7 
and Og as these mix under a scale evolution. 

Our current calculations use 1-loop matching 
factors for the 4-fermion operators. As explained 
in [Q , if the discretization of operators is such that 
Zp for the pseudoscalar bilinear is large, then the 



3. WILSON FERMIONS 

Our results with Wilson fermions are ex- 
ploratory. They have been obtained at just (3 — 
6.0, albeit on large lattices and with high statis- 
tics [||. The goal has been to understand sys- 
tematic errors, in particular the question of bad 
chiral behavior of matrix elements induced by the 
mixing with wrong chirality operators. 

The general form of the kaon matrix elements, 
as predicted by xPT, is 



K°(p f ) 0\K°(pi)) 



(8/3)/ 2 



+ Pm% +7i>» -Pf + 



/C,phys 



Sim K + 8 2 m 2 K p l ■ p } + 8 z {p l ■ p f ) 2 + 

where we ignore chiral logarithms and terms pro- 
portional to (m s — md) 2 ■ The former are diffi- 
cult to distinguish numerically from the terms we 
include, while the latter we expect to be small, 
especially for the range of quark masses studied. 

For Bk, chiral symmetry predicts that a, (3, 
and <5i are zero. With Wilson fermions, the mix- 
ing with wrong chirality operators generates these 
terms and, in addition, the allowed terms 7, 62, £3 
get contributions that have to be eliminated. 
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Operator 


Q 


= 6.0 


(3 = 6.2 


a = 


B 3 7 /2 


1/a 


0.989(05) 


0.823(16) 


0.62(3) 


B 3 7 /2 


it /a 


1.085(06) 


0.903(14) 


0.67(3) 


R 3/2 


1/a 


1.240(06) 


1.030(16) 


0.77(4) 


o3/2 
^8 


it /a 


1.288(06) 


1.076(17) 


0.81(4) 



Table 1 

3/2 3/2 

Results for ( NDR, 2 GeV) and Bg< ( NDR, 2 GeV), at the physical kaon mass, using smeared oper- 
ators. The last column gives the result of linear extrapolation in a 2 . 



With 1-loop improved operators these artifacts 
are 0(a 2 s ), but nevertheless overwhelm the signal. 
Our approach is to first remove a, (3, 8\ by study- 
ing the momentum dependence, and secondly, 
roughly estimate the artifacts in 7,(52,^3 using 
the fact that a, /?, Si have to be zero. Putting 
all these together gives, without extrapolation in 
a, 1 

Bk(NDR, 2 GeV) = 0.74 ± 0.04 ± 0.05 (5) 

where the second error is an estimate of the resid- 
ual contamination due to bad chiral behavior. 

Results for the electromagnetic penguins oper- 
ators, in the NDR scheme at fi = 2 GeV, are [B| 



B 3/2 = 



0.58 ± 0.02(stat) + ° Q ° 3 (pert) . 



(6) 
(7) 



B 3/2 = 0.81 ±0.03 (stat) (pert). 

The "perturbative error" reflects the dependence 
of the results on the choice of a s used in the 
matching of continuum and lattice operators, and 
is comparable to or larger than the statistical er- 
rors. A recent calculation of these by the APE 
collaboration using non-perturbative matching 
coefficients and the Csw = 1 clover action sug- 
gests that the errors in the 1-loop mixing coeffi- 
cients may be far more severe ||. 

The final quantities we consider are Bg = Bf 
and the related parameter B§ as defined in p|. 
The matrix elements we require are for bs mesons. 
The best we can do with present data is to give 
the result for ~ m c , 

Bf (NDR, 1/a) = 0.80 ± 0.01 (stat), (8) 
B£( NDR, 1/a) = 0.94 ±0.01 (stat). (9) 

A second limitation is that these results are at 
fi = 1/a = 2.33 GeV, the scale at j3 = 6.0, be- 



cause the two-loop anomalous dimension matrix 
needed to run to 2 GeV has not been calculated. 

4. CLOVER FERMIONS 

The analysis of data with clover fermions is 
preliminary. The same 170 lattices used in the 
study with Wilson fermions || are analyzed with 
tree-level tadpole improved clover action (Csw = 
1.4785). At this point we have a few qualitative 
statements about the data. 

The statistical fluctuations in the matrix ele- 
ments of S, V, A, V, T operators are much larger 
compared to those with Wilson fermions. Curi- 
ously, these fluctuations cancel in the five opera- 
tors Of . . . 0$ discussed above. 

The 1-loop mixing factors for O 3 ^ 2 and O 3 / 2 are 
too large and the calculation fails. The reason is 
the increase in Zp with Csw- 

The dominant artifacts in Bk, a and (3, are 
roughly a factor of five smaller compared to Wil- 
son fermions. As a result, the lattice value of Bk 
using the 1-loop improved operator improves from 
—0.30 with Wilson fermions to 0.50 with clover. 
This suggests that a large part of chiral violations 
is an 0(a) effect. 
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